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Abstract 

Using MacCormack-type methods, a new class of highly accurate compact MacCormack-type 
schemes is derived which does not require a tridiagonal matrix inversion to obtain the spatial 
derivatives. Two example schemes are shown, and results of these schemes for three linear and 
nonlinear CAA Benchmark Problems are presented. 

Introduction 

Computational aeroacoustics is concerned with the time-accurate solution of flow and acous- 
tic phenomena over long periods of time. To accomplish this goal, high-order finite-difference 
schemes and optimized schemes have been developed [e.g.. Refs. 1-6]. 

There are two main classes of high-accuracy finite-difference schemes: explicit schemes and 
compact schemes. Explicit schemes employ large computational stencils for accuracy, while 
compact schemes use smaller stencils by using the flux derivatives as independent variables at 
each grid point. While compact schemes are more accurate than the equivalent explicit scheme, 
solving for each flux derivative requires a scalar tridiagonal or pentadiagonal matrix inversion. 
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Recently, a new class of high-accuracy explicit MacCormack-type schemes have been intro- 
duced for computational aeroacoustics. 6 Applying this methodology to compact schemes has 
resulted in a new class of compact MacCormack-type schemes which use three-point stencils 
while retaining high accuracy. 

Mathematical Formulation 

Consider the 1-D equation: 


dU BF 
Bt Bx 


= 0 


( 1 ) 


A general compact spatial derivative with a three-point stencil can be written to calculate die 
spatial derivative: 


a D, _ ! + ( 1 - 2a)D, + a + j 




( 2 ) 


where Dj is the spatial derivative of F at point i . 

The standard Pade’ compact scheme 1 has the coefficients: 



( 3 ) 


This gives fourth-order spatial accuracy. Since the value of the derivative at a point depends on 
the value of its nearest neighbors, a scalar tridiagonal matrix must be solved to obtain the spatial 
derivative. 

Equation (2) can be written in matrix form as: 

[A]{D} = [F]{F} (4) 

where [A] and [B] are tridiagonal matrices. To solve Eq. (4), the matrix [A] must be inverted: 

{D} = [Af'[B]{F} (5) 

In a MacCormack-type scheme, the derivative operator is split into forward and backward 
operators such that: 




{D F ) + {P B } 
2 


( 6 ) 


The derivative operators can be written: 


{D F } = [A F ] _1 [5 F ]{F} (7) 

{ D B } = [aV[B B ]{F} (8) 
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Substituting into Eq. (6), 


{£>} = i([A F ]'V] + [A*fV]){F} 


(9) 


Multiplying through, we obtain: 


[. A F ][A B ]{D } = \{[A B }[B F ] + [A F ][F*])F (10) 


Notice that Eq. (10) is only valid if: 


[A F ][A B ] = [A B ][A F ] (11) 

Equation (1 1) is true for all tridiagonal matrices with constant diagonals. 

To avoid a tridiagonal matrix inversion, we use bidiagonal matrices for the forward and back- 
ward operators. Writing out explicitly, we have: 

( 12 ) 

> - ( * B + mS)Fi + 

( 1 - c F )D F i + c F D F i + 1 = (13) 

1 ■ ( * F + ™ F)F i + mFp i + 
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Using this, we find that at point i: 


[A F ][A S ]{Z>}|. = 

/( ! - c F )D , , + ( 1 - a s )( 1 - c")D, + /( 1 - a)D , , , 


Matching the coefficients in Eq. (14) with those on the left side of Eq. (2) gives: 


F B 
c -a 


b 1 ± J\ -4a 
a 2 



On the explicit side, we find: 


1 / , B.F.y-, 

-((a k )F t _ 2 

+ (( 1 - c F )k B + ( 1 - a B )k F - a B (k F + m F ))F f _ j + 

( c F m B - ( 1 - c F ){k B + m B ) - ( 1 - a B ){k F + m F ) + a B m F )F i 
+ ((1 - c F )m B - c F (k B + m B ) + ( 1 -a B )k F )F i + , 

+ (c F m B )F i + 2 ) 


Matching the coefficients in Eq. (16) with Eq. (2) gives: 



1 F B n 

k = -m = 0 


k B - -m = -2|3 = -1 (17) 

Performing a Taylor series expansion of the forward and backward operators gives: 

D F = F, + A*^F„-(A*) 3 ^F^ + 0(A*) 4 

D S = F I -AAr^F„ + (4x) 3 ^F ro „ x + 0(Ax) 4 (18) 

Clearly, a fourth-order accurate central difference is recovered when the forward and backward 
operators are substituted into Eq. (6). Using the terminology of Ref. [6], this method is a 4/2 
scheme. This terminology refers to the order of the underlying fourth-order central difference and 
the leading error term in the first-order biased differences. 

Notice that, since bidiagonal matrices are used to calculate the derivatives, the local value of D 
may be found using the value of D on one side only. By sweeping in the proper direction, the val- 
ues of D may be found explicitly. 

By using an additional explicit point and following the method used in Ref. [6], the order of 
the forward and backward operators may be increased to third order, resulting in a 4/4 scheme. 
This scheme is defined as: 


F B 
c = a = 


1 

3 


III 


6 



( 19 ) 


Consequently, one extra point is used in the explicit operator to achieve the increase in order. The 
Taylor series expansion of the forward and backward operators is: 


D f = F,-(A*) 3 ^F„„ + 0(A») 4 

D a = F x + (A + O(Ax)' (20) 


Substituting the coefficients of Eq. (19) into Eq. (16) reveals that the 4/4 scheme reduces to 
the following fourth-order compact scheme: 



-l 



+ !»i + i D ei 

If Af 

36 F, + 2 9 F, + 1 





( 21 ) 


This illustrates one very important difference between the explicit MacCormack-type schemes 
given in Ref. [6] and the compact MacCormack-type schemes described here. For explicit 
schemes, the order of the forward and backward operators may be changed without changing the 
underlying central difference (i.e., the sum of the forward and backward operators for the explicit 
4/2 and 4/4 schemes is the same fourth-order central difference, giving identical dispersion per- 
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formance), while the compact 4/2 and 4/4 schemes have completely different dispersion perfor- 
mance. 

Notice that using a larger stencil on the forward and backward operators results in a larger 
equivalent stencil for the underlying difference. As shown in Ref. [8], this method can be used to 
derive a class of prefactored small-stencil high-order compact schemes. 

Time Marching Methods 

The original MacCormack scheme used a two-stage explicit time marching method: 


C/ (I) = U n -Atb F [F(U n )] 

U n + l m Uu n + U {l) -AtS B [F(U W )]) (22) 


This scheme is second-order accurate in time, and is fourth-order accurate in space using the 4/2 
compact-implicit scheme given in Eq. (15) for linear problems. To obtain fourth-order accuracy 
in space for nonlinear problems, one needs to permute the forward and backward differences (see 
Ref. [5]). 

In addition, three multistage Runge-Kutta-type time marching methods for MacCormack-type 
schemes will be described. These methods were derived for MacCormack-type schemes in [6], 
and will be briefly described here. Using Eq. (1), a general six-stage Runge-Kutta MacCormack- 
type scheme can be defined as 


fc (1) = -Atb F [F(U n )] 


] i r 
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h m = + a 2 /» (1) )] 

h (3) = -Ar5 f [F(£/" + a 3 /i (2) )] 
h (4) = -A/8 B [F(t/" + a 4 /i (3) )] 


( 23 ) 


/j (5) = -Af5 F [F(J7 n + a 5 /i (4) )] 


Table 1: Time Marching Schemes 



RK2 

RK4 

LDDRK4- 
6 Step 1 

LDDRK4- 
6 Step 2 

«2 

1 

1/2 

1/2 

0.353323 

<*3 

0 

1/2 

1/2 

0.999597 

a 4 

0 

1 

1 

0.152188 

«5 

0 

0 

0 

0.534216 

«6 

0 

0 

0 

0.603907 

Pi 

1/2 

1/6 

1/6 

0.0467621 

P2 

1/2 

1/3 

1/3 

0.137286 

P3 

0 

1/3 

1/3 

0.170975 

P4 1 

0 

1/6 

1/6 

0.197572 

P 5 

0 

0 

0 

0.282263 

P6 

0 

0 

0 

0.165142 

Cl 

1 

1 

1 

1 

C2 

1/2 

1/2 

1/2 

1/2 

C3 

0 

1/6 

1/6 

1/6 

c 4 

0 

1/24 

1/24 

1/24 

c 5 

0 

0 

0 

0.0162098 

c 6 

0 

0 

0 

0.00286365 
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/i (6) = -Af5 5 [F(t/" + a 6 /i (5) )] 


u n+l = tf%M (1) + M (2) +M (3) 

+ m ( 4 ) +m ( 5 > +m (6) 

F B 

where 8 refers to a forward spatial difference and 8 to a backward spatial difference. Nor- 
mally, the order of forward and backward differences are interchanged every step to avoid numer- 
ical biasing. 

The coefficients for the three time marching schemes are given in Table 1. The RK2 scheme is 

second-order accurate, while the RK4 and Hu, et. al’s 4-6 LDDRK methods 9 are both fourth- 
order accurate in a linear sense. However, for nonlinear problems, the RK4 scheme is fourth- 
order accurate while the 4-6 LDDRK method presented here reduces to third order accuracy. 
Notice that the 4-6 LDDRK method uses a two-step marching cycle; one step has four stages and 
the second step has six. 


Accuracy Analysis 

Consider a general MacCormack-type method, which is Jfcth order accurate in time. The 
underlying central difference is /th order accurate in space, while the biased differences are mth 
order accurate in space. An accuracy analysis using Eq. (1) will give: 

U(t + At ) = U(t) + ^ + V(At k MAx m+l ,Ax l ) (24) 

ot 

From Eq. (18), the 4/2 scheme will have a third-order time-space error term, while the 4/4 


”1 i n 
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scheme is a true fourth-order method. 


A wavenumber analysis of these schemes gives a more complete picture of their linear perfor- 
mance and stability. The numerical wavenumber for a general three-point compact derivative (Eq. 
(2)) is defined as: 


©Ax = 


, /f -i0Ax ,, , x , /0 Axv 

-i(ke - {k + m) + me ) 

, -iQAx - / . \ . i0A*v 

(ae + l-(a + c) + ce ) 


(25: 


where 


—K < 0Ax < n (26) 


For the standard compact scheme, the numerical wavenumber is given by: 


0Ac = 


sin (6 Ax) 

(I + 5 C °s(0Ac)) 


(27) 


Notice that the right side of Eq. (27) is purely real. This is a property of central-difference 
stencils. 

For the 4/2 MacCormack-type scheme, the numerical wavenumber for the forward biased dif- 
ference is: 


li 



e F A* 


sin (6 Ax) 


—p{— 1 + cos(OAx))^ 

V 3 


2 , 1 

U3 + 3 C ° S 


(0A*)) Q + |cos(0Ax)) 


(28) 


Similarly, the numerical wavenumber for the backward biased difference is: 


0 B Ax = 


-p(- 1 + cosCOAx))^ 
sin (0 Ax) V3 

^0 + ^ cos (0Ax)j 0 + |cos(0Ax)J 


(29) 


There are two interesting things to observe in Eqs. (28) and (29). First, the real (dispersive) 
part of both equations are equal and identical to the real part of Eq. (27). Second, the imaginary 



Wavenumber 


Figure 1: Dispersion Characteristics of Compact MacCormack-type Schemes 
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Dispersion error 


(dissipative) parts are equal and opposite. This is a property of MacCormack-type schemes. 


Likewise, the numerical wavenumbers for the 4/4 scheme are: 


-F 

e ajc = 


v 

5 2 

' ~l8 + 9 


sin (0 sin ( 20 A*)) 

(§ + ! cos ( eA *>) 

(0A*) + ^cos(20Aar)y 

(J + | COS ( 0A *)) 


(30) 



Figure 2: Dispersion Error Magnitude Comparison 
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Figures 1 and 2 show the dispersion performance of these schemes, and compares them to 
explicit schemes such as 4th and 6th order central differences and Tam and Webb’s 7-point DRP 

scheme. 3 Notice the improvement in the dispersion performance throughout the wavenumber 
range from the 4/4 scheme; unlike all explicit MacCormack-type schemes, it outperforms the 
scheme that it is based on. 



Points per Wavelength 

Figure 3: Dissipation Magnitude for Compact MacCormack-type Schemes 
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Figure 3 shows the dissipation performance of the two MacCormack-type schemes. Again, 
the reduction in dissipation in the resolved wavenumber range from the 4/4 scheme is evident. 

Performance and Stability of Compact MacCormack-tvpe Schemes 

The numerical wavenumber of the forward and backward operators of a MacCormack-type 
scheme may be written as: 


0 F = 0 C -/5 (32) 

e* = e c + ;5 

Notice that the dispersion relation 0 C is the same for both the forward and backward differ- 
ences; however, the dissipative term 5 is equal and opposite. 

Putting these definitions into the time-marching schemes described above, we obtain: 


U(t + At) = ( 1 + cj (-iAf)0 C + (33) 

c 2 (-iAf) 2 ((0 C ) 2 + 8 2 ) + 
c 3 (-jA/) 3 0 C ((0 C ) 2 + S 2 ) + 
c 4 (-t'A/) 4 ((0 C ) 2 + 5 2 ) 2 + 

-C 2 o 2 

c 5 (-iAr) 0 ((0 ) +8 2 ) + 
c 6 (-iA0 6 ((e c ) 2 + 5 2 ) 3 )f/(0 
= (G(Q))U(t) 


Using this equation, we can find the frequency response of the time-marching scheme, as well 
as its stability. For this stability analysis, the magnitude of G(0) must be equal to or less than one 
for 0 < 0 < n. In this work, 0 -> 0 and 0 = n, which measure the long and short wavelength stabil- 
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ity, are investigated analytically. As a matter of observation, if the scheme is stable for small 0 
and for 0 = it, it is almost always stable for all 0. However, there is no mathematical proof that 
this has to be true. • 

Table 2 shows the results for the stability analysis, both analytical and numerical. Interest- 
ingly enough, the results for several of the schemes show that the most unstable wavenumber is 
not at 0 -> 0 or 0 = 7t. In particular, the LDDRK 4-6 scheme had its most unstable wavenumber at 
approximately 0 = 0.9 l7t. 

Boundary Stencils for Compact MacCormack-tvpe Schemes 

H Effect of Boundary Stencils on Interior S cheme Performance 

The numerical performance of the boundary stencil for a compact scheme has a much larger 
effect on the stability and accuracy of the scheme than the boundary stencil for the equivalent 


Table 2; Stability Limits 


Scheme 

Time 

marching 

CFL limit 
(0 -> 0) 

CFL limit 
(0 -> n) 

CFL limit 
(numerical) 

4/2 

RK2 

1 

73 

1 

73 

0.577 


RK4 

(stable) 

1 

1.0 


LDDRK 

4-6 

(stable) 

.934684 

0.891 

4/4 

RK2 

j 

(unstable) 

1 

2 

(unstable) 


RK4 

(stable) 

73 

2 

0.851 


LDDRK 

4-6 

(stable) 

.80946 

0.747 


























compact scheme. 10,11 The reason for this is that the error from the boundary stencil can propa- 
gate many points into the computational domain. 

Taking the 4/2 scheme as an example, let us assume that we have an error £q from the bound- 
ary at the beginning of the backward sweep. Defining Eq as: 


E 0 — (D ) interior ~ (D )boundary 


(34) 


where the subscript ‘interior’ refers to the spatial derivative that the interior scheme would have 
obtained and the subscript ‘boundary’ refers to the derivative obtained by the boundary stencil. It 
can be seen that we are defining the error with respect to the spatial derivative that the interior 
scheme would have obtained rather than the exact analytical derivative. 



Points from Boundary 


Figure 4: Effect of Boundary Stencil Error on Interior Solution 
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Substituting into Eq. (12), we find that the error propagates inward from the boundary. The 


error for a derivative i grid points away from the boundary is: 


( D ) interior i — (D )boundary\i + 


( B V 


(35) 


or 


e , = 


f B V 

a 

i B i 

1 -a ) 


(36) 


One point of interest is that the error due to the boundary stencil used at the start of the sweep 
has a much greater effect on the computed derivative than that of the boundary stencil used at the 
end of the sweep. Figure 4 illustrates the effect of boundary stencil error on the interior deriva- 
tives for the 4/2 and 4/4 schemes. It should be noted that the 4/2 scheme is less affected by 
boundary stencil error than the 4/4 scheme. 

2} Boundary Stencil Description 

One-sided, five-point, explicit boundary stencils were used for the two schemes described. 
These stencils were obtained by matching the Taylor series terms to the fourth order for the for- 
ward and backward stencils of each scheme. For space reasons, the stencils are not given here; 
however, the derivation is straightforward. 

This is not the only option for boundary stencil specification. Another possibility is to define 


1 I 
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ghost points outside the computational domain and extrapolate data to these artificial points using 
fourth-order accurate extrapolation. In this way, the interior scheme can be used on the bound- 
aries. For the fourth-order accurate compact schemes given here, only one ghost point is needed. 

^ F.ffect of Multistage Time Stepping on Accuracy at Boundari es 

The boundary condition specification for MacCormack-type schemes is not straightforward, 
due to the dissipative error terms in the one-sided differences. For example, the standard 
symmetry plane condition is set by using the method of images to set ghost values outside the 
domain. This, however, is only true for the first stage of a MacCormack-type method. For the 
second stage, the method of images is not correct; instead, the ghost values must be set in such a 
way that the end result is symmetric. For a two-stage RK2 method, this may be possible; for a 
four- or six-stage method, the errors can mount quickly. As noted above, the errors at the 
boundary will instantaneously propagate far into the domain at each stage. Thus, correct 
boundary condition specification is critical to the success and stability of compact MacCormack- 
type schemes. 

Benchmark Test Prob lems and Results 

Three test problems were chosen to investigate the numerical performance of the two compact 
MacCormack-type schemes. These problems are from the First and Second CAA Work- 
shops. 12,13 

n 1 -D Linear Wave P ropagation 

The first problem tested is the first linear problem given in the first CAA Workshop. The 
problem asks for the solution of the 1-D linear convection equation at time = 400: 
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0 


(37) 


U'+U x = 

U(x, 0) = iexp^-ln(2)^j j 


where: 

Ax = 1.0 (38) 

-20 <x£ 450 

Boundary condition specification is straightforward in this problem. At x = -20, which is the 
inflow boundary, the derivative of U is set to zero. At x = 450, which is the outflow boundary, the 
derivative of U is calculated explicitly from the interior. 

Figures 5-7 show the results for the 4/2 scheme used with each time-stepping method at vary- 
ing time steps. As noted in Ref. [6], the RK2 time-stepping method shows a tendency to affect the 



Figure 7: Performance of (4-6)-4/2 Scheme on Benchmark Problem 1 
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dispersion as the time step increases, while the RK4 method tends to be dissipative at higher time 
steps. The LDDRK 4-6 method also dissipates at higher time steps when combined with the 4/2 
scheme; however, this is more due to the scheme itself than the time stepping method. 

Figure 8 shows the results for the 4/4 method at a time step of CFL = 0.5 for both the RK4 and 
LDDRK 4-6 time stepping methods. The improvement in dissipation and dispersion from the 4/4 
scheme is evident, though the difference between the two time marching methods is small at this 
low time step. As the time step is increased, the LDDRK 4-6 method will show an improvement 
over the RK4 method. 

2) 1-D Shoc k Tube Problem 

The second one-dimensional problem solves the Category 2 shock tube problem from the 
First CAA Workshop. 12 The equations solved are the nonlinear Euler equations, written in con- 


0.40 h 


=> 0.20 


4k 

i\ 

I 
I 
I 

i 


... . 


J 4 


o Exact 

- 4-4/4 

- (4-6)-4/4 


- 0.20 


350.0 


400.0 

x 


450.0 


Figure 8: Performance of 4/4 Scheme on Benchmark Problem 1 




servative form: 


p 

d 

pu 

pu 

+ dx 

2 ► 
pu +p 

E 


. u(E + p ) . 


( 39 ) 


where 


P = <y-1)(£-jP« 2 ) 


(40) 


The equations are solved on the domain -235 < x < 235, with Ax = 1.0. The initial conditions 


are: 


u{x , 0) = 0 


P(x, 0) = 


4.4 x < —2 

2.7 + 1 .7 cos( (x+ /^ ) -2<x<2 
1.0 x>2 


p(x,0) = (YP) 1/Y 


(41) 


The boundary conditions used are Thompson’s 1-D characteristic formulation. 14,15 For 1-D, 


three characteristics are defined: 
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( 42 ) 


dp — du 
c * ■ d,- pc d, 


.2d 


dp 


e »- r a*»-af 


c 3 


= + pc— 

dt K 


where the overbar terms are mean values that are set as the initial values at the boundary. In this 
formulation, outgoing characteristics are computed using the interior scheme, while incoming 
characteristics are set to zero. For the inflow boundary (x = -235), c l and c 2 are incoming, and are 
set to zero. For the outflow boundary (x = 235), c 3 is incoming and is set to zero. 

Figures 9 and 10 show the results for density at time = 60 for this problem. In these figures, 
the expansion fan is shown on the left, with the contact surface in the middle and the shock on the 
ri gh t. Only part of the computational domain is shown in order to illustrate the performance of 
the different methods. 

In Fig. 9, the three time stepping methods are compared using the 4/2 scheme. All three use a 
time step of CFL = 0.25. The ability of all three methods to accurately track the discontinuities, 
with the RK2 time stepping method having more oscillations in the solution than the other two 
methods. 

In Fig. 10, the 4/4 scheme is tested using the RK4 and LDDRK 4-6 time stepping methods. 
At this s mall time step, both methods produce comparable results. Notice that, while the 4/4 
method does have more oscillations in the solution compared to the 4/2 method, the contact sur- 
face is more sharply resolved. This is due to the decreased dissipation and increased dispersion 
accuracy of the 4/4 method. 


1 I! 
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Figure 9: Effect of Time Stepping Method on 4/2 Scheme Solution of Nonlinear Benchmark 

Problem 2 (CFL = 0.25) 



Figure 10: Effect of Time Stepping Method on 4/4 Scheme Solution of Nonlinear Benchmark 

Problem 2 (CFL = 0.25) 
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The thir d test problem is from the Second CAA Workshop. 13 In this problem, a 2-D cylinder 
of radius R = 0.5 is located at the origin. At time = 0, an initial pressure pulse is specified: 


2 . 2 n 


P(x,y, 0) = exp 




The problem asks for the unsteady pressure data for 6 < t < 10 at three stations: 
(0,5),(5/V2 ,5/72), and (-5,0). 

The equations to be solved are the linearized Euler equations in polar coordinates: 


° Exact 

2-4/2 

— 4-4/2 

(4-6)-4/2 

a ° 4-4/4 

*• — -* (4-6)-4/4 


6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 

time 


Figure 11: Solution of Benchmark Problem 3 at (0,5). 
(CFL = 0.5) 
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( 44 ) 



The computational domain extends radially from R = 0.5 to R = 10.5, and azimuthally from 
0=0 to 0 = 7t. Three boundary conditions are used: a wall condition on the cylinder, an acoustic 
radiation condition in the far field, and a periodic condition at the azimuthal boundaries. 

The wall condition is based on Tam’s wall condition.^ This condition requires that the time 
rate of change of the normal velocity at the wall is zero: 

— r = j£ = 0 (45) 

dt dr 



Figure 12: Solution of Benchmark Problem 3 at (-5,0) 
(CFL = 0.5) 
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This condition is imposed by setting the normal derivative of pressure at the wall to zero for 


each sweep, while the other radial derivatives are computed normally. 
The acoustic radiation boundary condition is given as: 


L 

dt 


V r 

d 

v r 

1 

v r 

v e 

dr 

0 

r + T"i 

2r 

0 - 

* P < 


. p . 


. p . 


= 0 


(46) 


The radial derivatives at the exterior boundary are computed using one-sided explicit stencils. 

The periodic boundary condition requires some extra work. To start the sweep for the azi- 
muthal derivative, the one-sided boundary stencil is used. At the end of the sweep, a corrected 
value for the derivative is known, and Eq. (34) is used to update the interior points as required. 

For this calculation, a uniformly-spaced grid of 201 radial points and 301 azimuthal points is 

used, with a time step - CFL = 0.5. Since the grid does not have a point at (5/^2 ,5/ Jl ), data is 
only taken at the points (0,5) and (-5,0). 

Figures 11 and 12 show the results for these locations, compared to the exact solution. The 
results of both methods are good for this problem. At the initial transient in Fig. 1 1, it can be seen 
that the 4/4 scheme is less dissipative, but has the wave arriving slightly early. The 4/2 scheme is 
more dissipative, with a lagging dispersion error. At this relatively small time step, all time 
marching methods perform well. 


28 


IT 'IF 


Conclusions 


A new series of compact MacCormack-type schemes are presented which do not require tridi- 
agonal matrix inversions. Boundary stencils and boundary conditions are described, and results 
are shown for linear and nonlinear CAA benchmark problems. Stability and wave propagation 
performance are shown for the two schemes, using various time-stepping methods. The compact 
methodology gives a distinct performance advantage over previous explicit schemes, and can be 
easily added to existing MacCormack scheme solvers. 
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